Observation of a Triangular to Square Flux Lattice Phase Transition in YBa2Cu307 
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We have used the technique of smaU-angle neutron scattering to observe magnetic flux hues di- 
rectly in an YBa2Cu307 single crystal at fields higher than previously reported. For field directions 
close to perpendicular to the Cu02 planes, we find that the flux lattice structure changes smoothly 
from a distorted triangular co-ordination to nearly perfectly square as the magnetic induction ap- 
proaches 11 T. The orientation of the square flux lattice is as expected from recent d-wave theories, 
but is 45° from that recently observed in Lai.83Sro.i7Cu04+i . 

PACS numbers: 74.25.Qt, 61.12.Ex, 74.72. Bk, 74.20.Rp 



The technique of small-angle neutron scattering 
(SANS) from flux lines has a long and honorable record 
in measuring the properties of flux lines in superconduc- 
tors. However, it continues to bring new dividends, es- 
pecially in unconventional superconductors, since impor- 
tant information about the nature of the superconduct- 
ing state is often revealed by the flux line lattice (FLL) 
structure, for example 0,0,0,3. The diffraction pat- 
tern not only reveals the co-ordination and perfection of 
the FLL, and its correlation with the crystal lattice, but 
also the absolute intensity may be used to determine the 
actual spatial variation of the magnetic field within the 
mixed state and the values of the coherence length and 
penetration depth P, 0, 0| . In the simplest approxima- 
tion, flux lines would order in a regular triangular FLL; 
however, anisotropy of the Fermi surface or of the su- 
perconducting order parameter can cause distortions of 
the triangular lattice or transitions to other structures. 
The simplest situation in a high-K material is anisotropy 
of the magnetic penetration depth associated with effec- 
tive mass anisotropy 0; for example, the anisotropy in 
the ab plane of YBa2Cu307_5 (YBCO) leads at low val- 
ues of field to a corresponding distortion of triangular 
FLLs 01 • At the lower values of k in borocarbides, "non- 
local" effects are expected and observed 0, 0| to give 
a variety of FLL distortions and transitions. If the su- 
perconducting order parameter has a different symmetry 
from that of the crystal, this can again be revealed via its 
effects on the FLL structure, for instance in the p-wave 
superconductor Sr9Ru04 ^3- In general in d-wave 
superconductors [ij, , there is expected to be a ten- 
dency towards a square FLL as the field is increased and 
the anisotropic flux line cores overlap. This may be the 
cause of the FLL phase transition recently observed in 
overdoped Lai.83Sro.i7Cu04+d- (LSCO) at the compara- 
tively low field of 0.4 T 0. According to the FLL 
nearest-neighbour directions should lie along the direc- 
tions of the nodes of the order parameter, which would 
be at 45° to the Cu-0 bonds in the superconducting lay- 
ers. This is not, however, the orientation of the square 



FLL observed in LSCO The orientation may instead 
be controlled by band structure effects Q, even if the 
symmetry of the FLL is controlled by rf-wave effects. It 
has been suggested that an peak effect in magnetisation 
measurements on overdoped YBCO may be a signature 
of a continuous triangular-to-square FLL transition in 
this material at high fields ,15,] . However, others have 
suggested that there is a glass transition in this region 
[ig . Only by direct measurements may such suggestions 
be tested and the correlation between FLL and crystal 
lattice (or superconducting order parameter) determined. 

Our experiments were performed on the SANS-I in- 
strument at SINQ, PSI, Switzerland. Cold neutrons (8 A 
to 14 A, with a FWHM wavelength spread of 10 %) were 
coUimatcd over distances from 4.5 m to 15 m, depending 
on the field and hence g-range required. The diffracted 
neutrons were registered on a 128 x 128 x 7.5 mm^ multi- 
detector, which was similarly adjustable in distance from 
the sample. The undiffracted main beam was intercepted 
by a cadmium beamstop. A magnetic field of up to 11 T 
applied approximately parallel to the neutron beam, was 
provided by a cryomagnet with a field uniformity of 0.2 % 
over a 1 cm sphere. A variable temperature insert con- 
taining He heat exchange gas allowed sample tempera- 
tures from 1.5 K to 300 K. The sample was a 40 mg low- 
twin-density high-purity single crystal of YBCO grown 
in a BaZrOs crucible (l^ and oxygenated close to O7 by 
high-pressure oxygen treatment in order to reduce pin- 
ning by oxygen vacancies in the Cu-0 chains [l^ . It was 
therefore overdoped and had a Tc of 86 K. It was initially 
mounted with its c-axis parallel to the field direction. 
In order to satisfy the Bragg condition for each diffrac- 
tion spot in turn and hence establish the FLL structure, 
the cryomagnet and sample together could be rotated or 
tilted to bring the FLL Bragg planes to the appropriate 
small angles (~ 1°) to the incident neutron beam. 

In Fig. 1 is shown the FLL diffraction pattern ob- 
tained at the low field of 1 T. The most obvious fea- 
ture of this pattern is its fourfold symmetry which re- 
flects the average fourfold symmetry of the twinned or- 
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FIG. 1: FLL diffraction pattern at 1 T. The figure shows the 
counts on the SANS multidetector at 4 K (minus backgrounds 
obtained above Tc) summed over a range of angles between 
the field direction and the neutron beam. Noise at the centre 
of the picture has been masked. The cryostat was rocked by 
±1° about horizontal and vertical axes, ensuring that all spots 
in the diffraction pattern from the sample are detected. The 
{110} directions, corresponding to twin plane directions are 
vertical and horizontal in this picture. In all cases, the FLL 
was formed by applying the field above Tc and cooling. 



thorhombic structure of our YBCO sample. However, 
the FLL structure itself has triangular coordination, and 
the symmetry of Fig. 1 arises from four orientations of 
distorted triangular FLLs, present in different domains 
in the sample, as was first observed by Keimere< al. [l^ . 
The diffraction spots arising from these four triangular 
lattices are represented in Fig. 2. It appears that the dis- 
tortion of the individual lattices arises mainly from the 
a/b anisotropy present in each orthorhombic domain in 
the crystal |2Cll21| . This interpretation was confirmed by 
measurements on an untwinned sample which show 
diffraction spots distributed around an ellipse aligned 
with the a and b axes. The ratio of the principal axes 
of the ellipse should represent the anisotropy of the Lon- 
don penetration depth for Bd < -B < i?c2 The 
value we observe for the anisotropy ratio, "fab, in our 
sample is 1.28(1), whereas many estimates of this quan- 
tity are rather larger However comparable values 
to ours were obtained by measurements on a separate 
untwinned sample using neutrons Q and muons and 
torque magnetometry [23|. Our results are also corrob- 
orated by recent surface-sensitive measurements using a 
novel atomic-beam magnetic- resonance technique . It 
seems likely that the precise value of jab depends on the 
degree of perfection of the Cu-0 chains along the b di- 
rection [2^. The orientation of the triangular FLLs has 
been ascribed to pinning of a pair of spots, and hence 
planes of flux lines, to the twin planes |20ll2ll |. However, 
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FIG. 2: Schematic of diffraction patterns from four distorted 
triangular FLLs that together account for the pattern in 
Fig. 1. In each orthorhombic domain in the twinned crystal, 
there are two orientations of FLL (eg. (a) and (b)), derived 
by taking a regular hexagonal pattern and distorting it by 
the a/b anisotropy [Jj . The more intense pair of spots in each 
pattern is aligned with one of the {110} directions. The cen- 
tre figure is the superposition of the four FLL domains (a), 
(b), (c) and (d), which gives rise to the pattern in Fig. 1. The 
angles between reciprocal lattice vectors, a and /3, are defined 
for use in Fig. 5 



results reported later in this Letter also support the exis- 
tence of a correlation between the nearest-neighbour FLL 
directions and the directions of zeroes of the d-wave order 
parameter. 

In Fig. 3, we show diffraction patterns taken at higher 
fields. The data taken at 7T show a distortion of the 
triangular FLLs so that some of the weaker spots are 
closer to the strong spots, and others have moved towards 
the diagonals. There is clearly another source of distor- 
tion than pure a/b anisotropy. Finally at 11 T, the FLL 
has become almost exactly square, with the weak corner 
spots now playing the role of second order {1,1} spots of 
a square FLL instead of first order spots from a distorted 
triangular FLL. In order to investigate this steady change 
in the FLL structure with field, we rotated the crystal 
about the vertical axis in Fig. 3, so that the field was 5° 
from the c-axis. This was done in order to break the de- 
generacy between those FLL structures giving strong ver- 
tical diffraction spots (Fig. 2(a) and (c)), and those giv- 
ing strong horizontal spots (Fig. 2(b) and (d)). Within 
anisotropic London theory, this small angle of rotation 
should make a negligible change to the FLL distortion. 
As shown in Fig. 4., we found that at high fields the FLL 
structures giving horizontal spots were suppressed and 
instead only the structures depicted in Fig. 2(a) and (c) 
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FIG. 3: FLL diffraction patterns, as in Fig. 1 (but with a 
logarithmic intensity scale and smoothed to make the weaker 
spots more clearly visible), (a) at B = 7T, (b) at B = 11 T, 
showing the change in position of the weaker spots as the field 
is increased 



were observed. The advantage of this arrangement is that 
the pairs of spots near the horizontal axis in Fig. 4 could 
be observed easily, without being overlaid by the strong 
ones on the axis. This allowed us to measure accurately 
the spot positions and hence the FLL distortion. Never- 
theless, at high field, this pair of spots overlaps, but by 
assuming that the spot size is independent of field, we 
may estimate the angle between them even when they 
overlap. Further measurements of spot positions allow 
us to give a complete description of the FLL distortions 
versus field in terms of the angles between the FLL re- 
ciprocal lattice vectors. The results of this analysis are 
shown in Fig. 5. It is clear that the low field structure 
progressively changes with increase of field, although in 
our available field range the FLL never exactly reaches a 
perfectly square shape. This may partly be because the 
phase transition is at the extreme of our available field 
range, but is also clearly a result of the orthorhombic 
structure of YBCO. As depicted in the inset to Fig. 5, 
the a/b anisotropy of each domain must, on symmetry 
grounds, distort a square lattice to a rectangular one. 
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FIG. 4: FLL diffraction patterns with the field rotated 5° from 
the c-axis, to give only two FLL domains; (a) at _B = 7T, (b) 
at B = 10.8 T 



causing a slight splitting of the "square" spots from a 
twinned crystal such as ours. 

We have further investigated 0] the temperature- 
dependence of the FLL distortion shown in fig. 5. We 
find that with increasing temperature, the FLL structure 
changes more towards triangular. Thus the boundary be- 
tween square and triangular phases must curve up in field 
as temperature is increased. We would expect this if the 
triangular to square transition is due to d-wave effects, as 
the nature of the pairing becomes less important as 
becomes comparable with the magnitude of the gap. The 
shape of the phase boundary is similar to that seen in 
an overdoped sample by macroscopic measurements |l6j , 
but not the same as that proposed in ref. Unlike 
LSCO 0, the orientation of the FLL that we observe 
is aligned as expected from d-wave theories |0, 0| . It 
may be argued that twin planes, which are present in 
LSCO and YBCO are controlling the FLL orientation. 
To rule this out, measurements were also taken with the 
field at an angle to both sets of twin planes in our sample 
and the shape and orientation of the FLL was essentially 
unchanged. One should also note that the predicted dif- 
ference in free energy between the two orientations of a 
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FIG. 5; Field- variation at 5 K of two of the angles between 
reciprocal lattice vectors, a and depicted in Fig. 2. (Errors 
are comparable with the marker size.) All data were obtained 
with B parallel to the crystal c-axis except the data for a at 
fields greater than 6T, which were taken with B at 5° to c in 
order to resolve this angle more clearly (see text). A regular 
hexagonal lattice would have a = [3 = 60°, and an exactly 
square one a = 90° and j3 = 45°. Also marked by horizontal 
lines are predictions of anisotropic London theory Q for a and 
P, using a basal plane anisotropy, 'yab = 1.28, and assuming 
that one pair of spots is tied to the {110} directions. In the 
inset is shown the orthorhombic distortion from an exactly 
square pattern (exaggerated for clarity), expected in the two 
orthorhombic domains present in our crystal. 



square FLL is much larger than that between any trian- 
gular and the lower energy square orientation [l^. We 
also note that a similar correlation between FLL orienta- 
tion and probable direction of d-wave nodes has recently 
also been observed in CeCoIns Further support of 

the d-wave origin of the triangular-square transition in 
YBCO is the value of the transition field, which is a sim- 
ilar order of magnitude to the predicted 0.15i?c2 'iS^- We 
do not believe that in a large-K material such as YBCO, 
the transition is due to nonlocal effects as has been pro- 
posed for borocarbides 0| . We note that our hulk obser- 
vations of FLL structure are not in complete agreement 
with surface measurements by STM techniques [2^ Eoj. 
It is possible that the FLL adopts a different structure 
near the surface, or that the surface has different doping 
from the bulk. It is clear that the phase diagram is 
strongly doping-dependent. 

In conclusion, using small-angle neutron scattering, we 
have directly observed a change from triangular to square 
co-ordination of the flux line lattice as a function of mag- 
netic field in fully oxygenated YBa2Cu307. This phase 
transition is most naturally interpreted as a consequence 



of the d-wave character of the order parameter, which is 
expected to be more prominent at high magnetic fields, 
where the flux line cores begin to overlap. The orienta- 
tion of the FLL is as expected from an isotropic d-wave 
theory 0,0, unlike that in LSCO 0. It seems very 
likely from our measurements that at high fields, the FLL 
orientation, particularly when the FLL becomes square, 
is caused by the correlation between FLL planes and the 
directions of zeroes of the d-wave order parameter. It 
would clearly be of great interest to repeat these investi- 
gations in an untwinned crystal, where the effects of any 
pinning by twin planes would be completely removed, 
and the effect of the a/6 anisotropy on the "square" FLL 
would be crystal-clear. It appears that further investiga- 
tion of these phenomena will allow stringent tests of the- 
ories of the order parameter in the mixed state of high-Tc 
materials as a function of angle of field and doping. 

This work was performed at the Swiss Spallation 
Neutron-Source SINQ, Paul Scherrer Institut (PSI), Vil- 
ligen, Switzerland. 
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